Sauter-Schwinger like tunneling in tilted Bose-Hubbard lattices in the Mott phase 
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We study the Mott phase of the Bose-Hubbard model on a tilted lattice. On the (Gutzwiller) 
mean-field level, the tilt has no effect - but quantum fluctuations entail particle-hole pair creation via 
tunneling. For small potential gradients (long- wavelength limit), we derive a quantitative analogy 
to the Sauter-Schwinger effect, i.e., electron-positron pair creation out of the vacuum by an electric 
field. For large tilts, we obtain resonant tunneling related to Bloch oscillations. 
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Introduction. There are many striking analogies be- 
tween apparently very different areas in physics, which 
help us to unify our understanding of nature - as R. Feyn- 
nian said: The same equations have the same solutions. 
In this Letter, we establish such a quantitative analogy 
between high-field science and ultra-cold atoms in optical 
lattices. The Sauter-Schwinger effect [iHsf predicts that 
an extremely strong electric field may create electron- 
positron pairs out of the QED vacuum. As an intuitive 
picture, one can envisage this effect as tunneling of an 
electron from the Dirac sea into the positive continuum, 
leaving behind a hole (i.e., positron), see Fig. [1] Unfortu- 
nately, this prediction has not been experimentally veri- 
fied yet since it is hard to generate sufficiently strong elec- 
tric fields in large enough space-time regions. As we shall 
demonstrate below, the same equations govern particle- 
hole pair creation in a Mott insulator via tunneling due to 
a small potential gradient. The Mott insulator state we 
are considering can be generated by bosonic atoms in an 
optical lattice [4|, |5| for example, which repel each other 
and thereby become effectively pinned to the lattice sites 
(zero mobility). The benefits of this quantitative anal- 
ogy are two-fold. On the one hand, we can apply our 
knowledge of the Sauter-Schwinger effect [l|, y] to atoms 
in optical lattices and understand them better in this 
way. On the other hand, atoms in optical lattices may 
provide (via this analogy) an experimental approach to 
the so far unobserved Sauter-Schwinger effect. 
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FIG. 1: Sketch of the analogy: a) Dirac sea for E — 0, h) 
Sauter-Schwinger tunneling for _B 7^ 0, c) Mott state with 
energy gap AS, d) tunneling in tilted lattice. 



Sauter-Schwinger effect. Let us start with the Dirac 
equation describing relativistic electrons (ft = c = 1) 



7^(i9^-^g^^)*-m* = 0, 



(1) 



with charge q and mass m, propagating in an electromag- 
netic field Afj,. An electric field E can be encoded in the 
scalar potential E = V<i> or be generated by the vector 
potential E = dtA. Using the standard representation of 
the Dirac matrices 7'^ via the Pauli matrices cr and split- 
ting the bi-spinor ^ into upper and lower components 
^ = (V'-). ,tp-), the "square" of Eq. ^ yields 

{idt - q^f - (i V + qAf - m^ V± = iqE-CT i'^ . (2) 

For sub-critical electric fields m^ ^ qE, we may neglect 
the spin-dependent term on r.h.s. and obtain the Klein- 
Fock-Gordon equation for both spinors ip± separately. If 
the electric field is purely time-dependent E ~ E{t)ex, 
the temporal gauge $ = and E = dtA is most con- 
venient: After a spatial Fourier transform, we find that 
each mode k corresponds to a harmonic oscillator 
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(3) 



with a time-dependent potential. Replacing t by x, 
the above equation is equivalent to a one-dimensional 
Schrodinger scattering problem. A solution which ini- 
tially behaves as tp'ki^) '^ exp{—iujkt} will finally evolve 
into a mixture of positive and negative frequencies 
V'fc^nO o^ "fc cxp{~iujkt} + Pk exp{-fiiufct}. The Bogoli- 
ubov coefficient /3fe is related to the reflection amplitude 
in one-dimensional scattering theory and yields the prob- 
ability for electron-positron pair creation. For slowly 
varying fields E{t) with m? 3> qE, it can be derived via 
the WKB approximation. Since Eq. ([3]) corresponds to 
scattering above the barrier, the turning points are not 
real but complex. Hence the reflection (i.e., pair creation) 
probability is exponentially suppressed [l| demonstrating 
that it is a non-perturbative effect 
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where E denotes the maximum electric fleld. 



Bose-Hubbard model. Wc consider the Hamiltonian 



with the creation and annihilation operators aj^ and Ui, at 
the lattice sites /x and u, respectively. The tunneling ma- 
trix Tfj^i, G {0,1} encodes the lattice structure and J de- 
notes the hopping rate. The number of tunneling neigh- 
bors at any given site fi yields the coordination number 
f. Finally, U is the on-site interaction 



z = E.T, 



^f 



and V^ denotes the on-site potential with fi^ = aj^a^. 
This potential V^ at lattice position r^i will be analogous 
to the electric field E via V^{t) ^ qr^ ■ E{t). 

We assume an average filling of one boson per site 
{fifj} = 1. In order to have well-defined initial and fi- 
nal states, we envisage the following sequence: Initially, 
we have J = V^ = where the ground state simply reads 
l^in) = riull)/! (Mott insulator). Then we switch on 
the hopping rate J adiabatically (i.e., slow compared to 
the energy gap of the Mott insulator) such that we stay 
in the ground state. Next we introduce the lattice po- 
tential V^(t) which will be analogous to the electric field 
and enables particle-hole creation via tunneling. Again, 
we should do this slowly in order to avoid dynamical ex- 
citations y. Finally, we reverse this process and slowly 
switch off first V^(i) and then J. Thus the final ground 
state again reads Yl |1)^- If the final quantum state does 
not have exactly one particle per site |\E'out) 7^ Du |l)u 
we have a signature of particle-hole pair creation. 

In order to calculate this effect, we proceed along the 
lines of [6| and consider the time evolution of the density 
operator p of the whole lattice and introduce Liouville 
super-operators C^ and £^^ via 



idtp = 



H,p 



= 7 E ^^"'P + E ^^ 



(6) 
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As the next step, we derive the reduced density matrices 
for one lattice site p^ = Tr^{/5} via averaging (tracing) 
over all other sites v ^ p. and similarly for two and more 
sites /5^y = Tr^^{/5} etc. Then we separate the correlated 
parts via p^J^y = (f°" -I- p^pi/ for two sites, as well as 

A ;:corr i ;:corr ;: i ;;corr ;: i ;:corr ^ i A A A r^^r, 

PiLv\ - P^^yx + Pf_,y Px + P^\ Pv + Pv\ PiL + PilPvP\ etc. 
From Eq. ([S]), we obtain the effective evolution equation 
for one lattice site p via tracing over the rest v ^ p 



idtPtj. 



h^M 



K^fi 



oUp'^r + P.M} 



'-fiP/i : 



(7) 



where £^^ 



Cfi^i, + Cufi- In order to employ a con- 
trolled analytic approach, we consider the limit of large 
coordination numbers Z 3> 1. In this limit, the corre- 
lations obey the following hierarchy [6|, [D]. While the 
one-site density matrix p^ = 0{Z^) has entries of or- 
der one Tr^{/5p} = 1, the two-site correlations are sup- 
pressed via p™'' = 0{1/Z), the three-site correlations 



via p™;^ = 0(1/^2) and so on. Thus, Eq. © repro- 
duces the Gutzwiller ansatz to lowest order in 1/Z, which 
reads p^ = |1)^(1K 0{1/Z) = p^ + 0{1/Z) in the Mott 
state. Inserting this result and neglecting terms of order 
0(1/^2), we get for two sites 
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(8) 



Formally, this is an evolution equation for an infinite di- 
mensional matrix ffi^" , but fortunately it is sufficient 
to consider four elements only. Introducing local par- 
ticle and hole operators p^ = |l)u(2| and /i^ = |0) (1|, 
we find that their correlation functions /^j), = (/i^/i^), 

I^^y = {hlM, III = (Plhu), and f^l = (pIpu), obey for 
V ^ p a, closed linear system of equations 



_J_ 

~'z 



E 



T„ 



{idt + Vy, 



Z~ 






JV2 



T„ 



22 



Vu 



V2fil 



Vu) C 



12 



12\ 






(9) 



together with the symmetry f}^ = if SI,)*- Apart from 
this trivial relation, wc find an effective particle hole sym- 
metry f^ = ff^ (to first order in 1/Z). 

The Analogy. Eqs. ([9|) provide a complete set of equa- 
tions whose solution yields the number of particle-hole 
pairs created by the lattice tilt. However, instead of solv- 
ing them directly, we make the following trick. It turns 
out that we may recover (i.e., factorize) Eqs. Q to first 
order in 1/Z if we assume the following effective linear 
equations of the operators h^ and p^ 



V,- 



u 



Pf^ = 






idt -Vu + - 
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- ^ 2^ ^pi' 






+ V2K 
V2pu 



,(10) 



and exploit the initial conditions {hjj^h^)Q = S^^i, and 

{hlPu)o = {pihu)o = {piPu)a = in the Mott state. 
Note that these are effective equations of motion - they 
are not obtained by idthfj_ = [hfj_,H] etc. Nevertheless, 
the two-point functions in Eqs. behave as if the oper- 
ators hfj, and p^ evolve according to (jlOp. In this sense, 
solving Eqs. ^TU\\ yields the correct physics. 

As the next step, we consider the (continuum) limit of 
large length scales and approximate the tunneling matrix 



^ E T^Mt) « [1 + eV^ + OiV')] h{t, r) 



(11) 



by a long- wavelength expansion, where ^ — 0{1/Z) is 
the stiffness. Introducing effective scalar fields (j)± {t, r) = 
h{t, r) ± p{t, r), we thus find the set of equations 



(idt ~V)cj)± = 



1 



3 T 2\/2') V-U 



where P = J[l + ^V^ + C( V)]. As we shah see later, a 
very small tilt of the lattice will only induce a significant 
probability for pair creation if we are close to the critical 
point: U ^ J{3 + 2^2) + e with < £ < C/. In this case, 
the derivative part ^V^ can be neglected in the term 
(3 — 2a/2)2? and the coupled equations (fT2|) simplify to 



{idt 



vfh 



[r 



'cff-cff 



^off 



0{V')]4>+.{13) 



Up to short-distance corrections ©(V*), we obtain the 
analogue of the Klein-Fock-Gordon equation @ with an 
effective light velocity c^^j = ^(3JC/ — J^)/2 [12] and an 
effective mass TTicff which is related to the energy gap 



A£ = ^/.P -6JU + U^ = 2mcffc^ff 



(14) 



of the Mott insulator (in analogy to the energy gap 2meC^ 
between the Dirac sea and the positive energy solutions). 
Even if we are not close to the critical point, the above 
equation (J13p should remain a good approximation since 
the ratio (3 - 2^2)7(3 + 2\/2) « 0.03 < 1 of the two 
numerical factors in (|12p happens to be quite small. 

Particle-hole pair creation. For potentials of the form 
V^i(t) = qr^ ■ E{t), it is convenient to apply a Fourier 
transform after shifting the potential V^ (t) into the phase 



hf_t{t)e 



ip^dt'v^if) 



Y,hk{t)i 



.ik-r^ 



(15) 



and analogously for p^{t). Then Eqs. (|TOl) become 

tdthk = +^[iJTk-U]hk + V2JTkPk, 
idtPk = -\^[iJTk-U]pk-V2JTkhk, (16) 
where the Tk are now time-dependent 

Z^r^rr. ^k-(r -r ) rr. -^ I dt' [V^it')-VAt')] 

— 2^ Tfce''^ ('^^ '■■')=T^^e " . (17) 

k 

Inserting V^(t) — qr ^ ■ E{t) ~ dt[qr^ ■ A{t)], we see that 
the time dependence of T^ can be understood as replacing 
fc by k + qA{t). This is also known as Peierls substitution 
and is completely analogous to the gauge transformation 
in electrodynamics discussed between Eqs. ([T|) and ^. 
It also underlies the well-known Bloch oscillations. 
The most general solution of ([T5|) can be written as 



hk{t) = fi{Mk + Ik{t)Bk, 
Pk{t) = 9t{t)Ak+gk{t)Bk, 



(18) 



with suitable functions f^(t) and g^{t) and operators 
Ak and Bk- In the initial stationary regime where 
J = V^ = 0, we have hk{t) — /iJJ^cxp{i[/t/2} and 
Pk{t) = pJJ' exp{— i[/t/2}. Thus, w.l.o.g. we may iden- 



(12) tify h}^ with Ak and pj^" with Bk- This implies f^{t) 



exp{iJ7t/2} and g^{t) ex exp{— i[/i/2} initially with 
.9in ^ -/"in ^ *-*. During the time evolution according 
to ([T5|) . however, f^{t) and g^{t) will mix and posi- 
tive and negative frequencies will not stay separated - 
as discussed after Eq. ([3]) . In the final state (where again 
J := V^ = 0), we have hk{t) = /i^"* exp{i[/t/2} and 
Pk{t) = p™* exp{— iC/t/2} once more. Consequently, ini- 
tially and finally the positive or negative frequency com- 
ponents of the field 0^(t) = hk{t) + Pk{t) in Eq. P^ 
yield the initial and final particle or hole operators, re- 
spectively. As a result, a mixing of positive and negative 
frequencies of the 0+ -field in Eq. (|13p as given by the 
Bogoliubov coefficients ak and /3k directly corresponds 
to a mixing of particle and hole operators 



pr 



ak&^ + PkK 



(19) 



In Fourier space, the initial conditions discussed after 
Eq. (Uni) read {p\pk)in = {p\hk)in = {hlpk)in = and 
{h\hk)m = 1- Thus, inserting dH]), we get 



{pipk)c 



\i3k\ 



(20) 



This determines the on-site particle/hole probability via 
(2|pj2),,, = {plp,)out - EklPkWN = (0|pjO)„„t 
according to Eq. ([7]) and the aforementioned particle- 
hole symmetry f^l = f"^. As expected from the exact 
analogy between Eqs. ^ and P^ . we infer the same 
exponential scaling as in the Sauter-Schwinger effect ^ 



\Pk 



ex exp < — TT 



(A£)V4 + c,\fe 



IVylCeff 



(21) 



for small k and small and slowly varying lattice tilts. 
Note that this expression is non-perturbative in J and 
U , we only exploited the 1/Z-expansion. 

Bloch oscillations. The exact analogy established 
above applies to small lattice tilts and large length scales. 
For large potential gradients, the long-wavelength expan- 
sion is not applicable anymore and the lattice structure 
becomes important. For simplicity, we assume a square 
lattice in the following, with the gradient pointing along 
a lattice axis, but the results can easily be generalized. 
For large tilts, one obtains resonance effects which can 
be understood by considering the Fourier transform of 
Eqs. (m with the time-dependent Tk from Eq. dTT]) 

{idt-U + 2,JTk)fl^ - -V2JTki.fl'+fr + l), 
^dtfi' = idJl^ = V2JTk{fl^ - fl') , (22) 



where fl^ = Ul^)* ■ For a constant gradient VF, the Tk 
obey an oscillatory time-dependence due to the periodic- 
ity in fc ^ fe + t Vy , which is the basis of the well-known 



Bloch oscillations [5|, |8[ . For simplicity, we first consider 
the limit of small J which facilitates a perturbative solu- 
tion of Eqs. (1221) ■ To lowest oder in J, we obtain resonant 
growth if the potential difference /SV between two neigh- 
bouring lattice sites equals the gap A£ = U -\- 0{J) 



Ik — 



iJt 
272"^ 



-iUt 



0{J^ 



fk — 



JH' 



0{J^).{23) 



This resonance has a width of ^/2J and corresponds to 
tunneling to the nearest neighbouring site. In higher or- 
ders in J, we also obtain resonant tunneling to next-to- 
nearest neighbours for 2AF = A£ and so on. In contrast 
to the non-perturbative result ()2ip . this process is more 
similar to electron-positron pair creation in the pertur- 
bative multi-photon regime, see, e.g., [2|. 

Beyond perturbation theory in J, we may employ Flo- 
quct analysis [9| to find the various resonance bands 
which yield an exponential growth of the solutions to 
Eqs. ([22]). The associated Floquet exponent for the first 
resonance at AV ~ A£ reads J/(\/2Ay) and for the 
second resonance 2Ay = AS, it is 3J2/(4^/2[AV^]2) etc. 

Conclusions. In summary, we demonstrated that 
electron-positron pair production by a strong (and slowly 
varying) electric field and particle-hole pair creation in a 
slightly tilted Mott insulating Bose-Hubbard lattice are 
governed by the same equations. This quantitative anal- 
ogy is sketched in the following table: 



Sauter-Schwinger effect 


Bose-Hubbard model 


electrons & positrons 


particles & holes 


Dirac sea 


Mott state 


mass of electron/positron 


energy gap AS 


electric field E 


lattice tilt V^ 


speed of light c 


velocity CeS 



This analogy allows us to apply the vast machin- 
ery developed for the Sauter-Schwinger effect [l|-l3| to 
condensed matter theory. For example, small time- 
dependent variations of the potential gradient may signif- 
icantly enhance the tunneling probability |3| even if the 
rate of change is much smaller than the energy gap A£. 
These effects should be observable with atoms in optical 
lattices [Jl where the experimental technique now even 
allows the in situ detection of single atoms [10| . Finally, 
it would be interesting to study whether this quantitative 
analogy can be extented to the onset of dielectric break- 
down in the fermionic Hubbard model, see, e.g., 11|. 
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